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We determine the dispersion relations of fermionic quasiparticles in ultradegenerate plas- 
mas by a complete evaluation of the on-shell hard-dense-loop-resummed one-loop fermion 
self energy for momenta of the order of the Fermi momentum and above. In the case of 
zero temperature, we calculate the nonanalytic terms in the vicinity of the Fermi surface 
beyond the known logarithmic approximation, which turn out to involve fractional higher 
powers in the energy variable. For nonzero temperature (but much smaller than the chemi- 
cal potential), we obtain the analogous expansion in closed form, which is then analytic but 
involves poly logarithms. These expansions are compared with a full numerical evaluation of 
the resulting group velocities and damping coefficients. 
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Unscreened magnetostatic interactions in a degenerate Fermi gas lead to dramatic changes of 
the fermionic dispersion relation in the vicinity of the Fermi surface. At strictly zero temperature, 
there is a logarithmic singularity in the inverse group velocity, which leads to a breakdown of the 
Fermi liquid picture. This effect has been discovered in the context of a nonrelativistic degenerate 
electron gas by Holstein, Norton, and Pincus Q over thirty years ago, who found that it gives rise 
to an anomalous TlnT -1 behavior of the low-temperature specific heat (see also 0, 0, 0] ) ■ 

In deconfined degenerate quark matter the same effect is caused by unscreened chromomagnetic 
fields, since the nonperturbative magnetic screening is parametrically of the order g 2 T and thus 
vanishes in the low-temperature limit Although chromomagnetic screening may arise in the 
form of a Meissner effect in a color superconducting phase, the appearance of logarithmic terms in 
the quark self-energy at (resummed) one-loop order is also of importance in the case of color super- 
conductivity, since it leads to a a significant reduction of the magnitude of the superconductivity 
^ \ gap in a weak-coupling analysis 0, 0] ■ 

In the normal phase of degenerate quark matter, non-Fermi-liquid behavior leads to anomalous 
specific heat which because of the greater number of gauge bosons and the stronger coupling is 
comparatively large. Moreover, as has been shown in Ref. 0], the relevant logarithms are stable 
in the sense that they do not exponentiate into power-law behavior when higher loop orders are 
included as was previously assumed jsljlldj- However, an actual numerical evaluation requires to 
go beyond the leading logarithmic order calculations performed in the condensed-matter context 
EL Q]- This has been recently achieved for the low-temperature specific heat, where the 
scale of the leading temperature logarithm as well as subleading fractional powers of temperature 
were determined in Ref. [ill . H3 | , in a calculation which circumvented a complete evaluation of the 
fermion propagator. 

Non-Fermi-liquids effects have also been shown recently to significantly enhance the neutrino 
emission rate from normal quark matter E3l- However, this calculation involves the fermionic 
dispersion relations which have so far been known only to leading logarithmic accuracy. 

In this paper we shall close this gap and present a complete evaluation of the on-shell hard-dense- 
loop-resummed one-loop fermion self energy for momenta of the order of the Fermi momentum and 
above, both at zero and at small temperatures. 



II. FERMION SELF ENERGY ON THE LIGHT CONE 



The fermion self energy is defined through 

S- 1 (P) = S^ 1 (P) + Z(P), (1) 

where <So(P) = — (P) _1 is the free fermion propagator, and P^ 1 = (p°,p) with p° = iuj n + fj,, 
uj n = (2n + 1)ttT in the imaginary time formalism, and P^ = (E,p) after analytic continuation to 
Minkowski space. Without loss of generality we shall assume that [i > 0. 

With the energy projection operators Ap = i(l±7o7*]5*) we decompose X(P) in the quasiparticle 
and antiquasiparticle self energy, 

£(P) = 70 A+S + (P) - 7o A p S_(P), (2) 

and 

7o£'- 1 = a; 1 a+ + a: 1 a- (3) 

so that A± x = -\p° T (|p| + S±)]. 

The one-loop fermion self energy is given by 

s(p) = -JCjtY, I ^rs (p-Q)Y^ w (Q), (4) 



where is the gauge boson propagator. Following [20] we introduce an intermediate scale q*, 
such that m < 5* C |U, and we divide the q integration into a soft part (q < q*) and a hard part 

s + = 4 s) +4 h) . (5) 

For the hard part we can use the free gluon propagator, whereas for the soft part we have to use 
a resummed gluon propagator, see below. 

The hard contribution to S + on the light cone is most easily computed in covariant Feynman 
gauge, with the result 

4 h) = ^- (6) 

with = g 2 C/// 2 /(47r 2 ). Here q* enters only as a correction proportional to q*/fi, so that we 
can send q* to zero. Correspondingly we expect that in the soft contribution we should be able to 
send q* to infinity without encountering divergences, as will indeed be the case, but only after all 
soft contributions are added together. 

The leading and next-todeading contributions to the soft part of the fermion self energy for 
momenta of the order of \i or larger are obtained by a onedoop diagram where the fermion propa- 
gator is a bare propagator, but the gauge boson propagator is dressed by so-called hard-dense-loop 
(HDL) Q, m llfil Il7| self energies, so that the transverse and longitudinal parts of the gluon 
propagator are given by 0] 
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We shall consider zero or small temperature T <C p, for which the mass parameter in the HDL 
propagator is given by 



47T 2 



m 



(9) 



which is the asymptotic mass of the transverse modes, related to the Debye screening mass by 



m 



= mjj/2. 

On the light cone one finds the gauge-independent expression [l9( 







dt 



dk [S(ko -k)- 5(k + k)] 



x^2 ±sgn(fc )-p-qk-q / — Pt(<?0,^ 



2vr 



+ (±sgn(A;o) + kp 



2?r 



PL(qo,q) 



ko + qoT \E\ - ie 
1 + w fc (go) - n^fcp - /x) 
&o + 90 =F \E\ - ie 



"~2 ( 2 -«/( fc o-M) 



(10) 



where k = p — q and i£ = ±p. The distribution functions are given by n^qo) = l/(e <?0//T — 1) and 
nf(ko — fi) = l/Je^ 0- ^' 7 + 1]. pt and are the spectral densities of tranverse and longitudinal 
gauge bosons, respectively, 



PT,i(qo, q) = 2ImA TiL (g + ie, q). 



(11) 



We may use q <C \E\, k because of q < q* and \E\ > p. Depending on the sign of E, we can 
drop the term 5(ko + k) or the term S(ko — k) in Eq. ()1U|) . since its contribution is suppressed with 
~ q/E compared to the remaining contribution. Then we find for the soft contribution to the real 
part of £+ 



ReS 



(s) 



8vr 2 



/l " roc clq^ 
dt / — [(i -t 2 )PT(qo,q) + PL{qo,q)] 
-1 U-oo 71 



xV l+ n b (q ) - n f (E - p - gt) _ l_ ^ _ _ ^ _ ^ 



qo - qt 



(12) 



This quantity vanishes for E = p by symmetric integration. After performing the go-integration 
we therefore have 



ReS 



(s) 



g 2 c 



4vr 2 



dq q A 



dt (nf(E — p — qt) — rif(—qt)) 



(13) 



x [(1 - t 2 )Re A T (qt, q) + Re A L {qt, q)] . 
For Im S + (which receives no hard contribution) we find in an analogous way 

ImS + = -^-/ J dqq 2 J <ft ((1 - t 2 )p T (qt, q) + p L (qt, q)) 
x [1 + n b (qt) -n f (E- p- qt)] . 

The antiquasiparticle self energy £_ is obtained by inserting negative values of E in the expres- 
sions for and including an overall factor (—1). With p > we can then replace nj(E — p — qt) 
by 1. 



(14) 
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III. EXPANSION FOR SMALL \E fi\ AND SMALL T 



In this section we will perform an expansion of X+ in the region 

T ~ \E — n\ <C gn <C 



(15) 



where non-Fermi-liquid effects dominate. We will use the expansion parameter a := T/m, and we 
define A := (E - y)/T. From JTSJ) we have a«l and A ~ 0(1). 

In the part with the transverse gluon propagator we substitute q = ma}^z and t = a 2 ^ 3 v/z. 
After expanding the integrand with respect to a we find for the transverse contribution 



ReS 



00 

+(T) 



g Cfma I a 



7T^ 



3 A - 1 



av / ate . „ 

*1 7 a 2/3 k | (l + e^)(l + e A -^ 



2w 2 z(t; 2 7r 2 - 4z 6 ) 2/3 16v 4 z 3 (3w 2 vr 2 - 4z 6 ) 4/3 



U 2 7T 2 + 4z 6 (V 2 7T 2 + 4z 6 ) 2 



(v 2 vr 2 +4z 6 ) s 



(16) 



The z-integrations are straightforward. In the ^-integrals we may send the integration limits to 
=too. Using the formulae 



dv 



e x - 1 



dv- 



(1 + e v )(l + e x ~ v ) 
e 



T(a + 1) Li Q+1 (-e- A ) -Li a+1 (-e A ) Va > 0, (17) 



° A 1 



log \v\ 



7sA + g- ^Li Q+ i(-e _A ) - Li a+ i(-e A ; 



(1 + e v )(l + e A 

we find, neglecting terms which are suppressed at least with (m/q*) 4 , 
ReEjL = -g 2 C f m 

Alou I ^Ql) +1E \-® (u a+1 (-e- x ) - Li Q+1 (-e A ) 
am 6 Tt I oa V 



a=0 



(18) 



+ 



127T 2 

2 l/3 a 5/3 



o=0 



9x/3tt 7 /3 

2 2/3 a 7/3 



-20 



+ 



27vr 6 



r(|) ( L i 5/ 3(-e- A )-Li 5/3 (-e A )) 
T(l) (u 7/3 (-e~ x )-U 7/3 (-e x ) 
A (A 2 + rr } + 0(a 3 ) 1. 



27V3tt 11 /3 
8(24 - vr 2 )a 3 loga 2 2 



(19) 



In the longitudinal part we substitute q = rax and t = au/x. In a similar way as for the 
transverse part we find 



ReS 



(8) 



(L) = -9 Cjm 



a\ ( 2m 2 
^2 lo S 



(vr 2 -4)a 3 loga 2 2 



96vr 2 



A(A + 7r ) + 0(a 



(20) 



Turning now to ImS + we notice that it vanishes at E = fi only in the case of T = 0. For finite 
temperature, however small, there is an IR divergent contribution in the transverse sector 



Im £ 



00 



-(r)i£=M 



ln m 
4vr Air 



(21) 
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where the infrared cutoff may be provided at finite temperature by the nonperturbative magnetic 
screening mass of QCD. In QED, where no magnetostatic screening is possible, a resummation of 
these singularities leads to nonexponential damping behavior |2ll |. 
After subtraction of the energy independent part we have 



ImS 



(s) 



ImS 



CO i 



\E=/i 



a 2 C f q * f 1 

Z—£ dqq 2 dt(n f (E-fi-qt)-n f (-qt)) 
87r JO J-l 

x [(l-t 2 )p T (qt,q)+p L (qt,q)] . 
Following the steps which led to Eq. (|16j) . we find for the transverse contribution 



(22) 



ImS 



-(T) 



Im S 



(s) I 
-(T)\E=n 



g C fina / am 



2tt 



dv 



z 2 v 



+ 



16v z 



am 

3^4 



,1/3 



a I m 



a 2 /3|„| 



dz 



1 



(l+e v )(l + e 



X-i 



V 2 TT 2 + 4z 6 (V 2 7T 2 + 4z 6 ) 2 



2/3 16u 5 (f 2 7T 2 - 12Z 6 ) 4/3 

° + ( W 2 vr 2 + 4^)3 a +•• 



Using the formula 



dv 



7- T7- r — r\v\ sgn(a) 

(1 + e v )(l + e x ~ v ) 1 

-r(a + 1) [Li Q+1 (-e" A ) + Li a+1 (-e A ) + 2 (l - 2~ a ) ((a + 1) 



Va > 



we find in a similar way as above 



Im S 



(s) 



-en 



Im S 



00 



g z Cfml — — — log cosh ( — 



12vr 



+(T)IB=M 

(|) [Li 5/3 (-e-) + Li 5/3 (-e A ) + 2 (l - 2~ 2 / 3 ) C (f ) 

,7/3 . . . 

Li 7/3 (-e- A ) + Li 7/3 (-e A ) +2(1- 2" 4 / 3 ) C (I) 



-20 TTTrrd 



+ 0(a 3 



For the longitudinal part we obtain 



,(b) 



a 2 A 2 
64^/2 



+ 0(a 3 



(23) 



(24) 



(25) 



(26) 



We remark that the determination of the coefficient of the 0(a 3 ) terms in S + would require 
resummation of IR enhanced contributions along the lines of Ref. 42], App. A. 

Putting the pieces together, and using the abbreviation e = E — fj,, we obtain for the real part 



ReS. 



log 



4^2' 



KTh(e/T) 



+ 1 



2^V3 (T , 
45^73 Im 



5/3 



-20 



2 2 / 3 V3 (T , /£\\ 7/3 6144 - 256vr 2 + 36vr 4 - 9vr 6 / T , r e 



1897r n /3 V m 



TV 



864vr 6 



— /4 I tt; 
m \1 



log ( ^ 



(27) 
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where 



A (A) 
/2(A) 



exp 



1 d 

1-1E + (Li a+1 (-e- A ) - Li Q+1 (-e A )) 

3/5 



a=Q 



r(|) Li 5/3 (-e- A )-Li 5/3 (- e ; 



/3(A) = |r(f)(Li 7/3 (-e- A )-Li 7/3 (-e A )) 



3/7 



/ 4 (A) = A(A^+<) 



2MV3 



(28) 

(29) 

(30) 
(31) 



We note that the dependence on q* indeed drops out in the sum of the transverse and longitudinal 
parts. 

In the zero temperature limit (|A| — > oo) we have /i(A) — > |A|. If the temperature is much higher 
than \E-fi\ (i.e. A -> 0) we have / X (A) -> c := f exp(l - 7jB ) = 2.397357... and / 2 ,3,4(A) -> 0. 
For |A| 3> Co or |AJ <C Co we may approximate /i(A) with max(co, |A|), which is qualitatively the 
result quoted in It should be noted, however, that the calculation of Ref. |y] only took into 
account transverse gauge bosons, and therefore the scale under the logarithm and its parametric 
dependence on the coupling was not correctly rendered. 

For the imaginary part we find 



Im £+ — Im S- 



g Cfm 



24nm gi \T 

2 2/3 



+ 3 



2 l/3 

457r 7 /3 y m 

7/3 



fl2 



5/3 



where 



5i (A) = 2 log cosh 



647! {I 91 (£)) +20 i£w(^ 4 (f)) 



-) 



52(A) 
53(A) 

54(A) 



Li 5/3 (-e~ A ) + Li 5/3 (-e A )+ 2 (l - 2 



-2/3 



c( 



|A|, 



-r(f)(Li 



7/3V-e A )+Li 7/3 (-e A ) 



+ 2 fl-2- 4 / 3 V(l 



3/5 



3/7 



(32) 



(33) 

(34) 
(35) 

(36) 



In the zero temperature limit we have 5i(A) — * |A|. If the temperature is much higher than \E — fi\ 
we have 5i(A) — ► 0. 

Explicitly, our T = result reads 



IT=0 



M 2 
2E 



g 2 C f m 



12ir 2 m 

5/3 



log 



Ay/2: 



m 



7T £ 



+ 1 



+ 



l\S\ 

247rm 



•45^73- U J (sgn(e) -^' 

2 2 /3^3 



64^/2^^ 20 1897r u /3 
6144 - 256vr 2 + 36vr 4 - 97r b / : \, > 

m 



7/3 



J-i (sgn(e) + VSi) 
m 



m 



864vr 6 
11/3 



log 



0.928 m\ i7rsgn(e) 



(37) 
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where the scale of the last logarithm was determined by resumming IR enhanced contributions 
0. 

Apart from the first logarithmic term, the leading ima ginary parts contributed by the transverse 
and longitudinal gauge bosons were known previousl y Il9l . I23I . |24{ ] . As our results show, the damping 
rate obtained by adding these two leading terms [13, [23j is actually incomplete beyond the leading 
term, because the subleading transverse term of order |e| 5//3 is larger than the leading contribution 
from Tii,. 



IV. NUMERICAL RESULTS AND DISCUSSION 

At large values of E — [x or at large negative values of E, where one obtains the self energy 
of the antiquasiparticles, the soft contribution Re S s can be shown to vanish. At T = 0, where 
the imaginary part does not contain an infrared divergent contribution, ImS + approaches the 
constant 1 



E- 



lim ImE+(£7) = -g 2 C f m X 0.040534 .. . (38) 



T=0 



The resulting damping constant 7 = — ImS + is also that of the antiquasiparticles, which are of 
course far from their nonexistent Fermi surface for fi > 0. 

At intermediate energies \E — fj,\ > m, both the real and imaginary parts of E+ are nontrivial 
functions that we have evaluated numerically. The results are shown in Fig. Q for the two cases 
T = and T = m. Since at finite temperature the imaginary part of S contains a (constant) 
infrared singular contribution, we plot j(E — (j,) — 7(0) instead. The latter function is even with 
respect to its argument, resulting in a cusp at E = fj, for T = 0, while at finite T the damping 
contribution vanishes quadratically at E = fi. The real part is an odd function with respect to 
E — jj,. Again, E = fj, is a nonanalytic point at T = 0, but analytic at finite T. 



1 The numerical constant in Eq. H38H agrees with the value given with two significant digits in Ref. |2^l (taking 
into account the different normalization), and it is very close to, but not in complete agreement with, the value 
(2k + w 0.040854 quoted in Ref. |l<|. 
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FIG. 2: (v 1 — l)/(g 2 Cf) as a function of log 10 ((-E — //)/m) at zero temperature. 



At T = 0, the consequence of the nonanalyticity at E = fj, is that the group velocity dE/dp 
together with the residue in the propagator vanishes. The group velocity is determined by 



1 + 



2E 2 



dE 



(39) 



The numerical result for v g 1 



1 is given in Fig. El together with the series expansion for small 
E — ju| following from the results of the previous section. As one can see, this expansion converges 
well only for \E — fj\ <C m, and it turns out that the logarithmic contribution 



v~ x {E- l i) = l + 



12vr 2 



In 



4^/2 



m 



TT E 



3 

Ml + 2 



+ 0| ,'£Z£) 
m / 



2/3 



(40) 



is already a rather good approximation up to the point where one should switch to the leading order 
result M^/(2E 2 ) = g 2 C f /(M 2 ) 0.013£ 2 C/ that is relevant for larger values of m < \E- fj,\ < ft. 

At small finite temperature the results of the previous section show that the growth of v~ l for 
E — > [i is limited by 



g 2 Cr . c' n m _ 



(41) 



with d = 4V2e 5 / 2 /(irc ) ~ 9.15016. In this case the group velocity at and above the Fermi surface 
can be approximated by 



l/r, \ 9 2 C f f / 9.15 m 

[E - fi) « 1 + ^2^2 max ) mm I ln — f — 



,ln 



3.07m 



E 



Ml 



3 
: 2 



(42) 



This result for the group velocity of the quasiparticle excitations is for example of direct rel- 
evance for the calculation of the neutrino emission from normal degenerate quark matter [25| 
which is enhanced by non- Fermi- liquid effects As was shown recently in Ref. the neu- 

trino emissitivity involves two powers of a s ln(m/T) from the quasiparticle group velocities, which 
over compensate the single power of a s ln(m/T) in the specific heat that counteracts in the cooling 
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rate. For a numerical evaluation one evidently needs to know the constants under these logarithms. 
Interestingly enough, the constants under the log's in (1421) are much larger than the constant un- 
der the log arising in the specific heat determined in Ref. [lH Il2| as log(0.282m/T) — for more 
discussion see Ref. 12211 . 



V. CONCLUSIONS 



In this article we have computed the fermion self energy in an ultradegenerate relativistic 
plasma. For small \E — fi\ and small T we have obtained a perturbative expansion of £+ beyond 
the leading logarithm that is responsible for non-Fermi-liquid behavior. We found that dynamical 
screening leads to fractional powers in this series, which are analogous to the fractional powers 
in the anomalous specific heat in normal degenerate quark matter 11, 3]. Furthermore we have 
performed a numerical computation of the self energy and the group velocity for larger values of 
\E — n\. Our results provide an important ingredient for quantitative calculations of non- Fermi- 
liquid effects such as the computation of the enhanced neutrino emissivity of ungapped quark 
matter Il3l. 
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